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LONG-TERM  GOALS 

My  goal  within  the  EuroSTRATAFORM  program  is  to  understand  the  creation  of  the  preserved 
stratigraphic  record  on  continental  shelves  and  slopes  as  the  product  of  physical  processes  acting  with 
spatial  and  temporal  heterogeneities.  I  have  been  using  numerical  models  to  provide  insight  into  the 
formation  and  preservation  of  stratigraphic  sequences  at  margins.  My  goal  has  been  to  obtain  a 
quantitative  understanding  of  the  interactions  of  environmental  parameters  and  their  influence  on 
stratal  architecture  and  facies  distribution.  I  wish  to  be  able  decipher  the  stratigraphy  on  margins  to 
read  the  geologic  record  of  the  past  and  predict  future  stratigraphy. 

OBJECTIVES 

EuroSTRATAFORM  is  an  opportunity  to  extend  stratigraphic  model  applicability  to  new  margins  and 
to  expand  the  capabilities  of  our  products.  The  Rhone  and  Adriatic  margins  are  similar  to  the  earlier 
field  areas  in  that  they  are  mid-latitude  clastic  margins,  but  the  stratigraphy  is  highly  affected  by  the 
3D  geometry  of  the  basins.  The  forcing  provided  by  these  boundary  conditions  provide  a  challenge, 
but  are  also  necessary  for  modeling  realistic  systems.  My  aim  is  to  quantitatively  determine  the  system 
response  of  margins  to  different  forcing  functions  sufficiently  to  be  able  to  both  predict  stratigraphy 
and  invert  observed  sequence  architecture  for  geologic  history. 

APPROACH 

I  have  used  numerical  models  as  a  tool  to  provide  insight  into  the  fonnation  and  preservation  of 
stratigraphic  sequences  at  continental  margins.  To  study  these  margins  I  will  apply  a  multi-pronged 
approach.  Together  with  Italian  and  French  colleagues,  I  have  been  modifying  and  apply  our  existing 
2D  stratigraphic  models  to  serial  sections  of  the  two  EuroSTRATAFORM  margins.  This  will  test  the 
portability  of  the  models.  It  will  also  aid  in  evaluating  the  relative  along-  and  across-strike  transport 
and  tectonics,  and  their  influences  on  sequence  architecture.  John  Swenson,  with  assistance  from 
Chris  Paola,  Juan  Fedele,  I  and  others  are  jointly  developing  a  3-D  sequence  stratigraphic  model. 
Together,  we  will  extend  our  time-averaged,  moving-boundary  model  for  continental-margin 
sedimentation  with  its  coupling  of  multiple  transport  regimes  from  2-D  to  3-D.  The  advantage  of  this 
modeling  approach  is  that  it  allows  for  a  systematic  exploration  of  the  margin’s  response  to  variations 
in  sea  level,  sediment  supply,  tectonic  subsidence,  and  wave  climate  over  longer  timescales.  I  am 
providing  flexural  component  and  parameterizations  for  the  field  areas.  I  will  also  make  use  of  3D 
flexural  modeling  and  backstripping  being  done  for  a  separate  project  on  the  Gulf  of  Lion  to  recover 
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estimates  of  tectonic  subsidence  and/or  initial  morphologies  of  paleosurfaces  that  are  necessary  for  the 
forward  modeling. 

WORK  COMPLETED 

The  2D  model  developed  during  STRATAFORM  (Steckler,  1999;  Steckler  et  ah,  2001)  has  been 
modified  in  order  to  apply  it  to  the  Adriatic  field  area.  The  primary  changes  have  been  to  include 
along-strike  transport  of  offshore  sediments  and  improvements  in  variations  of  the  sediment  flux  with 
sea  level.  The  2D  modeling  is  focusing  on  the  Gargano  Peninsula  region  where  a  succession  of  four 
depositional  sequences  (e.g.,  Fig.  1)  mainly  composed  of  regressive  deposits  have  been  mapped  by 
Ridente  and  Trincardi  (2002).  Together,  we  have  chosen  a  series  of  10  seismic  profiles  along  the 
Gargano  structure  for  modeling  (Fig.  1).  We  started  with  one  of  the  northernmost  profiles,  which  is 
relatively  unaffected  by  the  tectonics  of  the  peninsula  and  are  now  extending  the  modeling  to  the 
south. 


Figure  1.  Map  of  Gargano  Peninsula  region  with  bathymetric  contours.  The 
thickness  of  Sequence  1  is  indicated  by  shades  of  gray.  Heavy  lines  indicate  the 
positions  of  the  10  seismic  profiles  selected  for  modeling.  An  inset  map  shows  the 
position  of  the  Gargano  Peninsula  region  along  the  Adriatic  coast  of  Italy. 

Led  by  John  Swenson,  a  3D  stratigraphic  model  has  been  developed  (Swenson,  2003;  Swenson  et  ah, 
2003).  I  am  working  on  completing  a  3D  variable  rigidity  flexure  module  (Graindorge  et  ah,  2004)  to 
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be  incorporated  into  the  code.  The  model  has  been  adapted  to  perform  simplified  simulations  of  the 
Adriatic  Sea  using  parameterizations  that  I  provided. 

A  new  parameterization  of  fluvial  dynamics,  which  includes  channel  belt  and  overbank  deposition,  has 
been  developed  by  Juan  Fedele  and  Chris  Paola  of  the  University  of  Minnesota.  This  provides  an 
improved  estimation  of  the  sand/mud  ration  and  thus  the  mean  grain  size  for  non-marine  strata.  This 
model  has  been  incorporated  into  the  multigrain-size  version  of  Sequence. 

RESULTS 

A  model  has  been  developed  that  reproduces  the  main  features  of  the  sequences  north  of  the  Gargano 
Peninsula  using  line  COS  124  (Fig.  2;  Steckler  et  al.,  submitted).  This  model  is  being  used  as  a  basis 
for  modeling  the  sections  farther  south  that  have  been  affected  by  the  tectonics  around  the  Gargano 
Peninsula.  In  the  profiles,  sequences  are  composed  primarily  of  regressive  deposits  separated  by 
regional  unconformities  (Ridente  and  Trincardi,  2002).  The  geometry  of  the  units  is  such  that  the 
regressive  clinoforms  flatten  out  at  their  seaward  termination  into  relatively  planar  strata  (Fig.  2). 
Attempts  to  simulate  the  sequences  were  frustrated  by  the  lack  of  thinning  of  the  sequences  seaward  of 
the  clinoform  rollovers.  Instead,  the  clinoforms  flatten  in  the  latter  stages  of  each  sequence.  This 
contrasts  with  the  development  of  clinoforms  and  a  depositional  shelf  break  that  is  characteristic  result 
of  most  conceptual  and  numerical  models  of  sequences.  We  therefore  modified  the  sequence  model  to 
account  for  the  dramatic  changes  in  sediment  pathways  that  happen  with  sea  level  fluctuations  in  the 
shallow  northern  Adriatic  Sea.  During  times  with  high  sea  level,  such  as  the  present,  sediment  from 
the  Po  and  smaller  Apennine  Rivers  are  transported  southwards  along  the  coast  by  marine  coastal 
currents  and  storms  (Cattaneo  et  al.,  2003).  At  times  of  low  sea  level,  fluvial  transport  of  an  enlarged 
Po  River,  in  which  the  Apennine  Rivers  are  captured  as  tributaries,  discharges  directly  into  the  MAD 
basin.  This  produces  a  reciprocal  sediment  supply  pattern  where  the  coastal  supply  at  high  sea  level 
produces  prograding  clinoforms,  but  their  supply  cuts  off  as  the  exposure  of  the  northern  shelf  at  low 
sea  level  switches  supply  to  the  fluvial  system,  which  discharges  into  the  260-m  deep  MAD.  In  the 
model,  the  clinoforms  flatten  as  nearshore  sediment  supply  decreases  and  is  replaced  by  the  deposition 
in  the  MAD  during  sea  level  low  stands.  The  inclusion  of  these  changes  in  sediment  supply  along  the 
profile  as  a  function  of  sea  level  enabled  us  to  obtain  a  good  fit  to  the  overall  stratigraphic  architecture. 
Additional  sediment  pathway  effects  during  sea  level  cycles  affect  the  sections  farther  south  due  to  the 
presence  of  structural  highs  such  as  the  Tremiti  Ridge.  Another  consequence  of  the  sediment  rate 
changes  was  that  the  decreasing  flux  during  sea  level  falls  aids  erosion  of  the  previous  transgressive 
deposits.  Finally,  the  model  required  uplift  increasing  towards  land  similar  to  estimates  reported  in  the 
literature  (e.g.,  Mastronuzzi  and  Sanso,  2002).  I  have  also  assisted  John  Swenson  in  creating  a  three- 
dimensional  model  of  the  Adriatic  Sea.  In  the  model,  current  sweep  the  sediments  along  the  Italian 
coast  during  highstand  but  the  shoreline  crosses  the  Adriatic  perpendicularly  during  lowstand. 

A  problem  with  models  of  fluvial  deposition  has  been  integrating  the  different  dynamics  of  channels 
and  floodplains  for  a  spatially-averaged  cross-section.  The  new  model  of  fluvial  deposition  predicts  the 
avulsion  frequency  of  channel  belts.  In  river  systems,  preserved  channel  belts  are  primarily  sand, 
while  the  overbank  deposits  are  primarily  muddy.  Avulsions  shift  the  channel  belt,  cutting  into 
previously  deposited  overbank  deposits  with  the  result  that  the  higher  the  avulsion  frequency,  the 
greater  the  average  percentage  of  sand  that  preserved  in  the  resulting  strata.  The  avulsion  frequency, 
among  other  things,  depends  upon  the  relative  mean  deposition  rates  in  the  channel  and  overbank 
regions,  the  inventory  of  sand  and  mud  in  the  fluvial  system,  subsidence  rates,  and  sea  level  change. 
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The  numerical  stratigraphic  model  can  readily  track  all  of  these  factors  and  calculate  their  variation 
both  downstream  and  through  time.  An  additional  result  is  an  estimate  of  the  channel  density  in  the 
preserved  strata.  Figure  3  shows  an  example  of  the  fluvial  section  of  a  model  run.  The  most  prominent 
feature  is  the  prediction  of  a  sandy  basal  above  the  sequence  boundary,  a  feature  common  in  outcrop 
and  coreholes.  In  addition,  subtler  variations  in  grain  size  with  fining  down  dip  and  at  maximum  sea 
level  rise  are  predicted,  again  features  consistent  with  field  observation  indicating  the  promise  of  this 
new  model. 


Figure  2.  The  top  shows  a  seismic  section  from  the  central  Adriatic,  located  on  map 
in  center  right.  Horizontal  scale  is  in  kilometers,  vertical  in  meters.  The  sequence 
boundaries  separating  the  four  depositional  units  representing  the  last  four  100  Ky 
sea  level  cycles  are  marked  by  colored  lines.  Below  is  a  simulation  of  the  section  at 
the  same  scale.  It  is  colored  acording  the  the  age  of  the  deposits  (indicated  on  sea 
level  curve  in  lower  left).  Unconformities  are  shown  as  white  lines.  The  stratigraphy 
is  composed  of  four  regressive  sequences,  labeled  1-4,  capped  by  Holocene 
transgressive  deposits  labeled  H.  The  model  fits  the  main  features  of  the  section  and 
the  geometry  of  the  strata,  namely  the  flattening  of  the  clinoforms  into  parallel 
bedded  strata  at  the  seawards  end  of  the  profde.  A  critical  part  of  matching  the 
depositional  patterns  is  reciprocal  sedimentation  pattern  caused  by  the  sediment 
source  switching  from  the  coastal  system  to  the  MAD. 
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Figure  3.  Cross-section  of  coastal  plain  from  numerical  model  colored  to  show 
lithology.  Horizontal  scale  is  in  kilometers,  vertical  in  meters.  The  top  section 
shows  timelines  in  blue  and  unconformities  in  red.  The  bottom  shows  the  identical 
section  without  the  lines  to  more  clearly  show  the  lithological  pattern,  especially  the 
coarse  lag  an  top  of  the  sequence  boundary. 


IMPACT/APPLICATIONS 

Sediment  supply  is  commonly  assumed  to  increase  during  sea  level  falls,  when  a  larger  area  is 
subareally  exposed,  and  decrease  during  sea  level  rises,  as  surface  area  is  flooded  and  sediments 
trapped  in  estuaries.  Along  the  Adriatic  coast,  we  find  the  opposite  effect  occurs  due  changes  in 
sediment  pathways  associated  with  sea  level  fluctuations.  This  alternative  sediment  supply  pattern  may 
be  common  along  margins  other  partially  filled  foreland  basins  where  sediment  transport  often  has  a 
large  along-strike  component.  It  may  also  occur  at  margins  where  rivers  discharge  directly  into  canyon 
systems  during  lowstands  such  as  both  the  Rhone  and  Eel  River  margins.  The  new  fluvial  model  holds 
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the  promise  of  better  prediction  of  lithology  and  channel  parameters  in  non-marine  strata  such  as  those 
that  underlie  many  continental  shelves. 

RELATED  PROJECTS 

I  have  been  funded  by  NSF  to  undertake  a  3D  backstripping  reconstruction  of  the  Gulf  of  Lions 
continental  margin  (Steckler  et  ah,  2003).  During  the  Messinian  Salinity  Crisis  (5.96-5.32  Ma),  the 
Mediterranean  was  cut  off  from  the  world  ocean  and  desiccated.  As  a  result  the  Gulf  of  Lions 
continental  was  heavily  eroded  and  the  margin  morphology  reduced  to  a  heavily  incised  seaward 
dipping  ramp  (Guennoc  et  al.,  2000).  Following  the  flooding  at  the  end  of  the  Messinian,  the  margin 
rebuilt  to  its  present  shelf-break  morphology.  The  60-100  km  of  progradation  during  the  Plio- 
Pleistocene  provides  a  stratigraphic  record  of  the  establishment  and  regrowth  of  the  margin  (Loll  et  al., 
2003,  in  press;  Steckler  et  al.,  2004).  The  3D  reconstructions  will  provide  reconstructions  of  the 
morphology  of  the  margin  for  time  slices  of  interest  to  EuroSTRATAFORM.  The  project  also 
involves  Greg  Mountain  and  Bill  Ryan  at  Lamont,  and  Serge  Berne  of  Ifremer  Marina  Rabineau  of 
IUEM  and  Christian  Gorini  of  USTL  from  France.  A  new  variable-rigidity  3D  flexure  code  has  nearly 
been  completed  for  this  project  (Graindorge  et  al.,  2004).  It  will  be  adapted  for  use  in  the  3D 
stratigraphic  model. 

REFERENCES 

Cattaneo,  A.,  A.  Correggiari,  L.  Langone,  F.  Trincardi,  2003.  The  late-Holocene  Gargano  subaqueous 
delta,  Adriatic  shelf:  Sediment  pathways  and  supply  fluctuations,  Marine  Geology  193,  61-91 

Graindorge,  D.,  M.  Steckler  and  L.  Lavier,  2004.  3D  deflection  of  a  lithospheric  plate  with  a  variable 
flexural  rigidity  and  under  a  variable  sedimentary  load  for  3D  backstripping:  preliminary  tests 
Submitted  to  AGU  Fall  Meeting. 

Guennoc,  P.,  Gorini,  C.  and  Mauffret,  A.,  2000.  Histoire  geologique  du  Golfe  du  Lion  et  cartographic 
du  rift  oligo-aquitanien  et  de  la  surface  messinienne.  Geologie  de  la  France,  3:  67-97. 

Lofi,  J.,  C.  Gorini,  M,  Rabineau,  S.  Berne,  G.  Clauzon,  P.  De  Clarens,  T.  Dos  Reis,  G.S.  Mountain, 
W.B.F.  Ryan,  M.S.  Steckler.  C.  Fouchet,  Plio-Quatemary  prograding  clinoform  wedges  of  the 
Western  Gulf  of  Lions  continental  margin  (NW  Mediterranean)  after  the  Messinian  Salinity  Crisis, 
Marine  Geology,  198,  289-,  2003. 

Lofi,  J.,  C.  Gorini,  S.  Berne,  G.  Clauzon,  A.T.  Dos  Reis,  W.B.F.  Ryan,  M.S.  Steckler,  Erosional 
processes  and  paleo-environmental  changes  in  the  western  Gulf  of  Lion  (SW  France)  during  the 
Messinian  Salinity  Crisis,  Marine  Geology,  in  press,  2004. 

Mastronuzzi,  G.  and  Sanso,  P.  2002.  Holocene  uplift  rates  and  historical  rapid  sea-level  changes  at  the 
Gargano  promontory,  Italy.  J.  Quaternary  Sci.,  17,  593-606. 

Ridente,  D.  and  F.  Trincardi,  2002,  Eustatic  and  tectonic  control  on  deposition  and  lateral  variability  of 
Quaternary  regressive  sequences  in  the  Adriatic  basin  (Italy),  Marine  Geology  184,  273-293. 


6 


Steckler,  M.S.,  High  resolution  sequence  stratigraphic  modeling:  1.  The  interplay  of  sedimentation, 
erosion  and  subsidence,  in  J.  Harbaugh,  L.  Watney,  G.  Rankey,  R.  Slingerland,  R.  Goldstein  and  E. 
Franseen  (eds.)  Numerical  Experiments  in  Stratigraphy,  SEPM  Memoir  62,  1999,  139-149. 

Steckler,  M.,  D.  Graindorge,  J.  Baztan,  M.  Rabineau,  C.  Gorini,  J.-L.  Olivet,  D.  Aslanian,  W.B.F. 

Ryan  and  G.  Mountain,  2004.  2-D  backstripping  of  the  Golfe  du  Lion  passive  margin  from  shelf  to 
basin:  insights  on  the  subsidence  and  paleo-geome tries.  Submitted  to  AGU  Fall  Meeting. 

Steckler,  M.S.,  J.  Lofi,  G.S.  Mountain,  W.B.F.  Ryan,  S.  Berne,  and  C.  Gorini,  Reconstruction  of  the 
Gulf  of  Lions  margin  during  the  Messinian  Salinity  Crisis  EGS-AGU-EUG  Spring  Meeting,  Nice, 
France. 

Steckler,  M.S.,  G.  Parker,  P.  Wiberg,  D.  Swift,  J.  Swenson,  C.  Reed,  L.  Pratson,  S. Parsons,  C.  Paola, 

A.  Niedoroda,  J.  Locat,  H.  Lee,  M.  Garcia,  S.  Fan,  J.  Carey,  2001.  Sequence4:  An  integrated 
stratigraphic  model  for  continental  margins,  at  “Formation  of  Sedimentary  Strata  on  Continental 
Margins”,  AGU  Chapman  Conference,  Puerto  Rico,  June  17-19. 

Swenson,  J.B.,  Modeling  clinoform  development  as  a  three-dimensional,  moving-boundary  problem, 
ComDelta:  Comparing  Mediterranean  and  Black  Sea  prodeltas,  26-28  October  2003,  Aix-en- 
Provence,  France 

Swenson,  J.B.,  A.B.  Murray,  C.  Paola,  M.S.  Steckler,  A  three-dimensional,  moving-boundary  model  of 
shelf  clinofonns,  EGS-AGU-EUG  Spring  Meeting,  Nice,  France,  2003. 

PUBLICATIONS 

Mountain.  G.,  R.L.  Burger,  D.L.  Orange,  L.F.  Pratson,  D.  Monteverde,  D.  Goldberg,  H.  Delius,  C.S. 
Fulthorpe, ,  C.M.  McHugh,  M.S.  Steckler,  K.G.  Miller,  Chapter  10:  The  long-term  record  on 
continental  margins,  in  C.  Nittrouer,  J.  Austin,  M.  Field,  M.  Steckler,  J.  Syvitski  and  P.  Wiberg 
(eds.),  Continental-Margin  Sedimentation:  Transport  to  Sequence,  125pp.,  submitted,  2004. 

Pratson,  L.,  J.  Swenson,  A.  Kettner,  J.  Fedele,  G.  Postma,  A.  Niedoroda,  C.  Friedrichs,  J.  Syvitski,  C. 
Paola,  M.  Steckler,  E.  Hutton,  C.  Reed,  M.  Van  Dijk  and  H.  Das,  Modeling  Continental  Shelf 
Formation  in  the  Adriatic  Sea  and  Elsewhere,  Oceanography,  submitted,  2004. 

Pratson.  L.,  P.  Wiberg,  M.  Steckler,  D.  Cacchione,  J.  Karson,  E.  Mullenbach,  J.  Swenson,  C.  Nittrouer, 
B.  Murray,  G.  Spinelli,  C.  Fulthorpe,  D.  O’Grady,  G.  Parker,  N.  Driscoll,  R.  Burger,  C.  Paola,  D. 
Orange,  M.  Wolinsky,  M.  Field,  C.  Friedrichs  and  J  Fildelez,  Chapter  8:  Seascape  evolution  on 
clastic  continental  shelves  and  slopes,  in  C.  Nittrouer,  J.  Austin,  M.  Field,  M.  Steckler,  J.  Syvitski 
and  P.  Wiberg  (eds.),  Continental-Margin  Sedimentation:  Transport  to  Sequence,  170pp., 
submitted,  2003. 

Steckler,  M.S.,  D.  Ridente,  F.  Trincardi,  2003.  Modeling  the  Stratigraphy  and  Tectonics  Around  the 
Gargano  Peninsula  on  the  Adriatic  Coast  of  Italy.  COMDELTA  Conference,  Aix-en-Provence, 

Oct.  26-28,  2003. 

Steckler,  M.S.,  D.  Ridente,  F.  Trincardi,  Modeling  the  Stratigraphy  and  Tectonics  Around  the  Gargano 
Peninsula  on  the  Adriatic  Coast  of  Italy,  AGU  Fall  Meeting,  2003, 


7 


Steckler,  M.S.,  D.  Ridente  and  F.  Trincardi,  2004.  Control  of  sequence  geometry  north  of  Gargano 
Peninsula  by  changing  sediment  pathways  in  the  Adriatic  Sea,  Continental  Shelf  Research, 
submitted. 

James  P.M.,  L.  Pratson,  P.  Wiberg,  M.  Steckler,  M.  Garcia,  R.  Geyer,  C.  Harris,  E.  Hutton,  J.  Imran, 
H.  Lee,  M.  Morehead,  G.  Parker,  Chapter  11:  Prediction  of  margin  stratigraphy,  in  C.  Nittrouer,  J. 
Austin,  M.  Field,  M.  Steckler,  J.  Syvitski  and  P.  Wiberg  (eds.),  Continental-Margin  Sedimentation: 
Transport  to  Sequence,  197  pp.,  submitted,  2003. 


8 


